Spin dynamics of current driven single magnetic adatoms and molecules 



O 

(N 



F. Delgado, and J. Fernandez-Rossier 
Departamento de Fisica Aplicada, Universidad de Alicante, San Vicente del Raspeig, 03690 Spam 

(Dated: June 30, 2010) 

A scanning tunneling microscope can probe the inelastic spin excitations of a single magnetic 
atom in a surface via spin-flip assisted tunneling in which transport electrons exchange spin and 
energy with the atomic spin. If the inelastic transport time, defined as the average time elapsed 
between two inelastic spin flip events, is shorter than the atom spin relaxation time, the STM 
current can drive the spin out of equilibrium. Here we model this process using rate equations and 
a model Hamiltonian that describes successfully spin flip assisted tunneling experiments, including 
a single Mn atom, a Mn dimer and Fe Phthalocyanine molecules. When the STM current is not 
spin polarized, the non-equilibrium spin dynamics of the magnetic atom results in non-monotonic 
dl/dV curves. In the case of spin polarized STM current, the spin orientation of the magnetic atom 
can be controlled parallel or anti-parallel to the magnetic moment of the tip. Thus, spin polarized 
STM tips can be used both to probe and to control the magnetic moment of a single atom. 
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I. INTRODUCTION 

A single magnetic atom is arguably the smallest sys- 
tem where the spin can be used to store classical and / or 
quantum information. Therefore, there is great inter- 
est in probing and manipulating the spin state of a sin- 
gle atom or a single molecule in a solid state environ- 
ment. Examples of this are single Phosphorous donors in 
Silicon,^ nitrogen- vacancy centers in diamonds, 2 ~ 5 single 
Mn atoms in II- VI 6 -^ and III-V— semiconductors, and sin- 
gle magnetic adatoms in surfaces^—. Whereas in most 
cases the spin of the single atom is probed by optical 
means, the possibility of coupling the spin of as single 
atom to an electrical circuit is particularly appealing. 

Tremendous recent experimental progress has made it 
possible to probe the spin of a single and a few atoms 
deposited in conducting surfaces by means of scanning 
tunneling microscopes^ - — There are two complementary 
techniques that afford this: spin polarized STM and spin 
flip inelastic electron tunnel spectroscopy (IETS). The 
working principle of spin polarized STM is spin depen- 
dent magneto resistance^ similar to that of tunnel mag- 
neto resistance junction: tunneling between two spin po- 
larized conductors depends on the relative orientation of 
their magnetic moments. Control of the spin orientation 
of either the tip or the substrate affords spin contrast 
STM imaging, 1 ! 

In the case of spin flip IETS, electrons tunnel from 
the STM tip to the surface (or vice versa) , and exchange 
their spin with the atom, so that they produce a spin 
transition, whose energy is provided by the bias volt- 
age. Thus, a new conduction channel opens when the 
bias voltage is made larger than a given spin transition 
[see Fig. 1(c), 1(d)]. This results in a step in the con- 
ductance as a function of bias and permits to determine 
the energy of the spin excitations, and how they evolve 
as a function of an applied magnetic field.— r 12 i 15 i 16 i 18 i 22 
When the atom is weakly coupled to its environment, the 
spin is quantized ) 23 ' 24 the spin transitions have sharply 



defined energies which can be described with a single ion 
spin Hamiltonian whose parameters can be inferred from 
the experiments ^ 12 ' 16 This is the case of Mn, Co and Fe 
atorns ^ 11 : 16 as well as Fe and Co PhthalocyanineS ) 12 ' 15 : 18 
all of them deposited on a insulating monolayer on top of 
a metal. Remarkably, spin flip IETS does not require a 
spin polarized tip to extract information about the spin 
dynamics. 

Both IETS and spin polarized STM are based upon the 
fact that the spin state of the atom affects the transport 
electrons, yielding a spin-dependent conductance. There- 
fore, we must expect that the transport electrons do af- 
fect the spin of the atom. This is the main theme of this 
paper. In the case of spin flip IETS, there are two rele- 
vant time scales. On one side, the inelastic transport time 
or charge time T q , which is defined as the average time 
elapsed between two inelastic spin flip events . On the 
other side, the magnetic atom spin relaxation time, T±. 
In the T q » T\ regime, the transport electrons always 
interact with a atomic spin in equilibrium with the envi- 
ronment. As a result, the occupation of the spin states 
is bias independent and the conductance is expected to 
have flat plateaus in between the inelastic steps. In the 
T q « T\ regime this is no longer the case, the current 
drives the atomic spin out of equilibrium, so that the oc- 
cupations of the spin states are bias dependent. As we 
showed in a previous work^ for the case of a single Mn 
atom, this results in a modified conductance line-shape, 
with non-monotonic behaviour in between steps. In this 
work we give an extended account of these effects, and 
consider also the case of Mn dimers and FePc. 

Non-equilibrium effects become particularly appealing 
when either the tip or the substrate are also spin polar- 
ized. In this case we have current flow between two mag- 
netic objects, which is expected to result in spin transfer 
torque^ It has been proposed theoretically^ and inde- 
pendently verified experimentally— that the spin orien- 
tation of a single Mn atom can be controlled with a a 
spin polarized tip. In this paper we provide a thorough 
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FIG. 1: (Color online), (a) Scheme of the proposed setup: 
a magnetic STM tip and a magnetic adatom on a insulating 
monolayer deposited on a metal. A current flows through the 
adatom when a bias voltage is applied between tip and sur- 
face, (b) Energy level diagram of tip, atom and surface and 
typical microscopic process that gives rise to spin-flip tunnel- 
ing, (c) Equivalent circuit that accounts for the increase of 
the conductance schematically shown in (d). 



analysis of this effect and extend our study to include the 
effect of an external magnetic field, as in the case of the 
experiments^ 

The results presented in this paper are based on 
a phenomcnological spin-depedent tunneling Hamilto- 
nia n 22 i 25 ' 28 " — and are, in most instances, in agree- 
ment with existing experiments. In the case of a sin- 
gle magnetic atom with spin 1/2 it is possible to de- 
rive our spin-dependent Hamiltonian from a single or- 
bital Anderson model, by means of a Schrieffer- Wolff 
transformation ! 32 ' 33 Within this picture, the spin-flip as- 
sisted tunneling events would correspond to inelastic co- 
tunneling in the Anderson model. In the single orbital 
Anderson model the spin-flip channel can dominate the 
elastic channel, which can be even zero in the so called 
symmetric case. Further work 34 is in progress to gener- 
alize this picture to the higher spin case relevant to this 
paper. 

The spin dynamics of current driven nanomagnets in 
the Coulomb Blockade regime has been thoroughly stud- 
ied from the theory standpoin t 25 ' 35 " — . The systems 
studied include magnetic grain o 41 ' 45 , semiconductor— 
and Mn doped quantum dot a 36 ' 40 , molecular magnets 
and magnetic molecule a 35 ' 42 ' 43 ' 46 . In contrast to most 
of the theory work published to date, here we can 
model the current driven dynamics of a quantum spin 
whose Hamiltonian parameters are accurately known 
from experiments 1 ^, making it possible to compare suc- 
cessfully theory and experiment. 

The rest of the manuscript is organized as follows. In 
Sec. II we present the model Hamiltonian for the mag- 
netic atom(s), the transport electrons, and their cou- 



pling, which accounts both for spin-assisted tunneling 
and Korringa like atomic spin relaxation due to exchange 
coupling with the electrodes. The transition rates and 
non-equilibrium dynamics leading to the current are an- 
alyzed in Sec. Ill CI In Sec. Ill we present the results of 
current driven spin dynamics under the influence of non- 
magnetic tip in three cases: the single Mn adatom, the 
Mn dimer and the FePC molecule. In Sec. IV we discuss 
the case of a spin polarized tip and analyze in detail the 
case of a single Mn adatom. In Sec. V we present our 
main conclusions and discuss open questions. 

II. THEORY 

A. Hamiltonian 

In this section we present the phenomenological Hamil- 
tonian, its microscopic justification, the rate equation ap- 
proach for the atom spin dynamics, including both spin 
relaxation and spin driving terms, and the calculation 
of the current. The system of interest is shown in Fig. 
QJa). We use a model Hamiltonian which describes the 
system of interest split in 3 parts: tip, substrate and the 
magnetic atom(s}2 5 . 

H = H T +n S + Hgpin + V. (1) 

The first two terms describe the tip and surface surface: 
n T + U s = e ^( fc ) c iL, c fc<"j = (2) 

where cj[ creates and electron in electrode ?y = T, S, 
with momentum k and spin a defined along the spin 
quantization axis, n. Unless stated otherwise, we take 
ft parallel to the magnetization of the tip, which is a 
static vector in our theory. Since we consider a non- 
magnetic surface, we have e CTi s(fe) = eg(k). All the re- 
sults of this paper are trivially generalized to the case of 
a non-magnetic tip and a magnetic surface. 

The spin of the magnetic adatom(s) is (are) de- 
scribed with a single ion Hamiltonian, exchange cou- 
pled to other magnetic adatoms and to the transport 
electrons^.22,25 

Hspi„ = Y. [ D &® + E (X 2 « - S 'y^) 

i 

+ \ E JijS'aW-S'aV) + 9»b E S'(t).B. (3) 

The first term describes the single ion magneto- 
crystalline anisotropy, the second describes the inter- 
atomic exchange couplings and the third corresponds to 
the Zeeman splitting term under an applied magnetic 
field B. Here the prime denotes that the spin quantiza- 
tion axis is chosen with z' along the easy axis of the sys- 
tem, not along the magnetic moment of the tip, n. This 
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makes necessary to rotate Wspin when ft is not parallel to 
the easy axis. The value of the local spin S(i), the mag- 
netic anisotropy coefficients D and E, and the exchange 
coupling between atoms in the chain Jij, change from 
atom to atom and also depend on the substratei^ r 11 ' 22 ' 47 
In the following, we denote the eigenvalues and eigenvec- 
tors of "Hspin as Em and \M) respectively. 

We model the coupling of the magnetic chain 
with the reservoirs with the following Kondo-likc 
Hamiltonian* 2 ^* 2 ^ 2 ^ 



V= T x ^, a {i) T -^S a {i)c\ i(T c xl<T ,, (4) 

where i labels the magnetic atoms in the surface, A = 
(fc, rf) labels the single particle quantum numbers of the 
transport electrons (other than their spin a) , and the in- 
dex a runs over 4 values, a = x,y, z, and a = 0. We 
use and S a for the Pauli matrices and the spin op- 
erators in the n frame, while So = t^ ' is the identity 
matrix. T\ t \> ia for a = x, y, z is the exchange-tunneling 
interaction between the localized spin and the transport 
electrons, and potential scattering for a = 0. Attend- 
ing to the nature of the initial and final electrode, Eq. 
(U) describes four types of exchange interaction, two of 
which contribute to the current, the other two are crucial 
to account for the atom spin relaxation. 

B. Justification of the Hamiltonian 

The phenomenological spin models of Eqs. © and (0| 
capture most of the experimental results, as we show be- 
low. These models imply that that the magnetic atom is 
in a well defined charge state except for classically for- 
bidden fluctuations that enable tunneling from the tip to 
the surface. Fig. Hfb) shows a typical level alignment in 
which the spin model can be applied. The basic condi- 
tion is therefore, that the chemical potential of the elec- 
trodes must be far enough from the chemical potential of 
the central-quantized region. In this way, charge addition 
and charge removal are classically forbidden. Although it 
is outside the scope of this work, we claim that the quan- 
tum charge fluctuations that give rise to spin dependent 
tunneling are due to inelastic cotunncling. In the case of 
spin 1/2, the equivalence between the spin model, origi- 
nally proposed by Appelbaum 2 ^ and a single site Ander- 
son model was rigorously shown by Anderson 32 generaliz- 
ing the Schrieffer and Wolff transformation 3 - 3 - to the case 
of a single site coupled to two reservoirs. Within this pic- 
ture, the atomic spin is exchanged coupled to the trans- 
port electrons and the magnitude of the exchange is given 
by Tx^\', a =x,y, z — 5~ 1 V\V\>, where V kn is the hybridiza- 
tion between the Anderson site and the single particle 
state k in the electrode 77, and S is the energy difference 
between the Anderson level and the electrode Fermi en- 
ergy. This is the so called Kinetic exchange. Importantly, 



both electrode conserving and electrode non-conserving 
processes are included and their strengths are not inde- 
pendent, since they both depend on hopping matrix ele- 
ments V\ between the localized orbital in the atom and 
the extended orbitals in either the tip or the sample. In- 
terestingly, the Schrieffer and Wolff transformation 3 ^ also 
yields a spin-independent tunneling term which would 
yield the a = contribution in Eq. Q, and it corre- 
sponds to the elastic tunneling contribution. Within this 
Anderson-Kondo picture, the strength of the (a = 0) 
elastic channel and that of the spin dependent channel 
(a = x,y,z) are comparable and, in the so called sym- 
metric case, the elastic term vanishes identically. Thus, 
for spin 1/2 case, this picture can account for the large 
strength of the inelastic signal. The generalization to 
higher spin case, relevant for the experiments^ -12 ' 16 ' 27 
will be published elsewhere.™ 

Keeping these considerations in mind, and follow- 
ing Anderson^ we assume that Hamiltonian (Q} arises 
from kinetic exchange. The momentum dependence of 
Tx.\', a {i) can have important consequences in the con- 
ductance profile^ in an energy scale of eV, but it can be 
safely neglected in IETS. We thus parametrize 

T VlV >, a (i) = v v (i)v v -{i)%, (5) 

where vs(i) and vt(i) are dimensionless factors that scale 
as the surface-adatom and tip-adatom hopping integrals. 
Because kinetic exchange is spin rotational invariant we 
have % = T y = T z = |T|. Thus, Eq. (gj implies that 
the spin- assisted tunneling and the atomic spin relax- 
ations are both due to kinetic exchange, and Eq. ([S]) im- 
plies that their amplitudes depend on the tip-atom and 
surface-atom tunneling matrix elements. 



C. Rates and Master equation 

Our primary goal is to study transport and spin dy- 
namics. This is done considering V as a perturbation to 
the otherwise uncoupled magnetic atom and transport 
electrons. The quantum spin dynamics is described by 
means of a master equation for the diagonal elements of 
the density matrix, Pm, described in the basis of eigen- 
states \M) of %s P i n . The master equation is derived us- 
ing the standard system plus reservoir technique^, where 
the transport electrons act as a reservoir for the atomic 
spin(s). The master equation 49 reads 

^T = I2 Pm ' Wm '< m - Pm E Wm - m ' ' ( 6 ) 

M M' 

where Wm,M' are ^ ne transition rates between the atomic 
spin state M and M'. These rates can be written as 

Wm,m' = J2n,n' W m?m> > wner e W^ M , are the scattering 
rates from an atomic spin state M to M' in which a 
quasiparticle electron goes from electrode r\ to rj' as a 
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result of exchange process. They are given by: 

wv$.= v c: fe v M '/,(^)[i-^(^)] (7) 



M,M' ~ / j 
kk' .era' 



where f n (e) = 1/(1 + exp [—j3(e — n v )]) is the occupation 
probability in electrode r\ for electrons in equilibrium at 
chemical potential p^ and temperature T = l/(/cg/3). 

k'a-'M' * s ^ ne ra ^ e a ^ wri ich an electron in lead r/ with 
wavenumber k and spin a is scattered into a lead 77' with 
wavenumber fc' and spin a', with the impurity spin un- 
dergoing a transition between states M and M'. Quan- 
tum rates Vs are calculated at the lowest order in the 
electrode-chain coupling using Fermi Golden rule with 
the perturbation given by V (see appendix A for details) : 



■pn-tv 

1 kaMMa'M 1 
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2T a «„(i)v(i)^ lJ "'(i) 

a, z 

x<5 (e^(fc) + E M - e a , v ,{k') - E M ,) , (8) 
where we have defined the matrix elements 



S. 



MM' l a 



(i) = (M\S a (i)\M'). 



(9) 



The rates in Eq. ([7} describe 3 types of processes: 



1. Elastic processes, W M M , in which the state of 
atomic spin remains unchanged and a transport 
electron is transfered from one electrode to another. 
These processes are responsible for the elastic cur- 
rent and have no effect on the spin dynamics. The 
rates of the elastic processes scale with v^VgT§. 

2. Spin transitions W^ M ,. In these, a spin transition 
in the atomic spin is produced due to the creation 
or annihilation of an electron hole pair either in 
the tip or in the surface. These processes do not 
contribute to the current. At very small tempera- 
ture, the fastest process of this type is atomic spin 
relaxation: a spin transition from an excited state 
Em to a lower energy state Em' which results in the 
excitation of an electron hole pair in one of the elec- 
trodes. This spin relaxation process is very similar 
to the nuclear spin relaxation due to hyperfine cou- 
pling to conduction electrons in metals and to Mn 
spin relaxation in diluted magnetic semiconductors 
due to itinerant carriers^. At zero bias these pro- 
cesses dominate the atomic spin relaxation time T\ . 
The rates of the spin transition processes scale like 

4T 2 and 4T 2 . 

3. Spin flip assisted tunneling W 7 ^ M ,. In these pro- 
cesses, which contribute both to inelastic current 
and to the dynamics of the atomic spin, a trans- 
port electron goes from electrode rj to rf inducing 
a spin transition from state M to M 1 . The rates of 
the spin flip assisted tunneling processes scale like 



2-r2 

T U S> ■ 



The steady state solutions of Eq.© depend, in gen- 
eral, on the Hamiltonian parameters, the temperature 
and the bias voltage. We refer to the steady state so- 
lutions as Pm(V). At zero bias, the steady state solu- 
tions are those of thermal equilibrium. At finite bias, the 
Pm{V) can depart significantly from equilibrium depend- 
ing on the relative efficiency of the transport assisted spin 
excitations and the spin relaxation. Eq. © does not in- 
clude spin coherences. This approximation is good if the 
spin decoherence is faster than spin relaxation, which is 
known to be the case due to hyperfine coupling^ in Mn 
atom. However, future work should address this point 
more carefully, in particular when the magnetization of 
the tip n is not parallel to the single ion easy axis. 



D. Relevant parameters 

The behaviour of the system is characterized by the 
rates in Eq.([7]), which depend on a number of physical 
quantities like the temperature, the bias voltage V, den- 
sity of states at the Fermi Energy of tip and surface, pt 
and ps, the tip-atom vx(i) and surface-atom vs{i) hop- 
pings, the spin independent To and spin dependent T 
couplings. In this work we attempt to group the unknown 
parameters either in terms of dimensionless numbers or 
as experimentally accessible quantities. For that matter, 
we define the zero bias elastic conductance 



— v n \rr 1 2 2 

.9o = -rC-oPTPs \T \ X , 



(10) 



where Gq = 2e 2 /h is the quantum of conductance and 



22v T {i)vs(i) 



(11) 



is a parameter that quantifies the tip-surface transmis- 
sion through the magnetic atoms. The density of states 
at the Fermi energy for spin a in the electrode 77 are de- 
noted by p^ . We define the spin assisted conductance 
g s as 



9s 



C 2 9o 



where 



(12) 



(13) 



is the ratio of the spin-flip assisted and elastic tunnel 
matrix elements. We shall use spin polarization of the 
tip, defined as 



PTt — PT\ 



Pt\ + Pt\ 

Another important parameter is the ratio 

v T (i) 



r[i) 



vs(i) ' 



(14) 



(15) 
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which decreases as the tip is retracted from the surface. 
In most instances, we shall have r(i) < 1. As a gen- 
eral rule, the processes that drive the magnetic adatom 
out of equilibrium are proportional to v^v^ whereas the 
processes that cool the spin down (if k^T < \eV\) are pro- 
portional to + Vg. Thus, the non-equilibrium effects 
are higher as 72. = r 2 /(l + r 2 ) increases. Therefore, cur- 
rent will be increased without changing the applied bias 
V . In contrast, the inelastic ratio £ and the magnitude 
of the tip polarization Vt are not so easy to control. 



E. Current 

The calculation of the rates for a tunnel event in which 
a transport electron goes from one electrode to the other, 
inducing a spin transition between states M and M' 
(where M could be equal to M' in the elastic channel), 
permits to obtain an expression for the current in terms 
of the steady state solutions of the master equation 



MM- 



M,M' 



(16) 



where WjJfM' are the scattering rates from state M to 
M' induced by interaction with a quasiparticle which is 
initially in reservoir rj and ends up in r)', given in equation 
0. We adopt the convention that positive bias voltage 
V > means electrons flowing from tip to surface, see 
Fig. Efb). Thus, we have: 

eV = /is - Mt (17) 
with e the value of the electron charge with its sign. 

1. Current for non-magnetic tips 

The expression for the current in the case on non- 
magnetic tip and substrate can be written as the sum 
of two terms, elastic and inelastic, I = Iq + IlN given by 
the expressions^ 



= 9oV, 

where go is given by Eg. (1101) and 



IjN = 7T i-(Alf,M' 
Go — ' 

u MM', a 

with 



eV) 



iM,M' 
*a,TS 



(18) 



Pm (t0 (19) 



<w = iE^)MOWa(OK)- (20) 

x i 

Here we have introduced the current associated to a sin- 
gle channel with energy E±m,m' = Pa/ ~~ Pm'> an d bias 
V 

i±(A + eV) = — (G(A + eV) ± G(A - eV)) , (21) 
e 



with G(uj) = oj (l — e _,3LJ ) 1 . The curve i-, is odd in the 
bias, whereas i+, relevant in the case of magnetic tips dis- 
cussed below, is even. In contrast, di^/dV is even and 
di + /dV is odd. In this non-magnetic tip case, the elas- 
tic current provides no information about the spin state 
whatsoever. In contrast, the inelastic steps in conduc- 
tance arise from Eq. (|19j) and permit to extract informa- 
tion about the spin transition energies, Am,M' and spin 

matrix elements Sf^ s . The basic effects of the elastic 
and inelastic terms in the conductance can be understood 
in terms of an equivalent electric circuit schematically 
shown in FigOJc) an d (d). For low enough voltage, the 
only channels that can conduct current are the elastic 
ones, while the inelastic channels remains close. In this 
situation the (inelastic) switch is open. When the voltage 
is increased such as inelastic channels are open (switch 
is closed), these new channels contributes to the current, 
leading to a smaller resistance, see Fig.QJd). 



2. Current for magnetic tips 

In the case of a magnetic tip, the current has 3 con- 
tributions, I = Iq + I mr + 1 in- This result is different 
from the non-magnetic case on two counts. First, the 
elastic case has a magnetoresistive term, so that current 
is now proportional to the relative orientation of the av- 
erage adatom spin and the magnetic moment in the tip 
n: 

h + Imr = go [1 + 2C(S,,ts)V t ] V, (22) 

where 

(S z ,ts) = Ys P M{V)v T (i)vs{i){M\S a (i)\M) (23) 

M,i 

is the average magnetization along the z axis, that we 
take parallel to the magnetic moment of the tip. As we 
show both below and in Ref. H^, both Pm(V) and the 
average atom magnetization depend on voltage. Impor- 
tantly, the magnetoresistive contribution to the elastic 
current makes it possible to track changes in the single 
atom magnetization experimentally. 

The second difference with the non-magnetic tip arises 
in the inelastic current, which is now given by the 
expression 



2o.2<) 



lin 



M [i-(A M , M ,+eF)£ 

M,M' 



3 a,TS 



+P T i+(A 



V)S xy (M,M')\p M (V). (24) 



MM' 



The new term in the second line involves the matrix ele- 
ments 



3 xy (M,M') = 2Im 



qM,M' qM' ,M 
S x,TS a y,T,S 



(25) 



As opposed to the standard inelastic current [Eq. ([15])]. 
which gives rise to steps of equal height for positive and 
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negative bias, the Vt dependent term of the inelastic 
conductance, proportional to i + , yields steps at the exci- 
tation energies of opposite sign as the polarity of the bias 
is reversed. Both the elastic and inelastic term propor- 
tional to Vt can produce a dl/dV which is not an even 
function of bias. 



III. NON MAGNETIC TIP 

In this section analyze the implications of a non- 
equilibrium population distribution when a finite bias 
is applied between tip and surface. These effects will 
be more relevant when the current through the sys- 
tem increases (by increasing the coupling to the elec- 
trodes). Next, we will study these effects in three differ- 
ent systems: the Mn monomer, Sec. lIII Al the Mn dimer, 
Sec. IIIIBl both deposited on a CU2N surface and the iron 
Phthalocianine molecule, FePc, deposited on an oxidized 
Cu(lllO) surface, Sec. IIICl 



A. Mn monomer 

Let us consider first the case of a single Mn adatom in 
CU2N, which has been widely studied experimentally 9 -^ 
and theoreticall y 22 ! 29 '^^— The spin of the Mn atom 
in this environment is S = 5/2. The parameters of the 
single ion spin Hamiltonian have been determined exper- 
imentally to be D = -0.039 meV, E = 0.007 meV^ and 
g = 1.98. 9 Since E « \D\ we can limit our qualitative 
discussion to the case E — 0, so that the eigenstates of 
%Spin are also eigenstates of S z (numerical simulations 
will be done with E = 0.007 meV and do not change 
qualitatively). In the absence of applied magnetic field 
and at temperatures much smaller than the zero field 
splitting 4|D|, the equilibrium distribution is such that 
the two ground states, S z = ±5/2, are equally likely and 
the average magnetization is zero. At an energy of 4|D| 
above the ground state level, we find a couple of degen- 
erate excited states, with S z = ±3/2. Finally, the two 
states with S z = ±1/2 are found at 6|D|. 

From the experiments, performed at low current^ the 
experimental dl/dV lineshape is piecewise constant with 
two steps symmetrically located eV — ±4Z?. This is ac- 
counted for by the equilibrium theory^ 2 - As we show in 
Fig. [2] and also in our previous work^ non-equilibrium 
effects modify the dl/dV lineshape. In particular, the 
dl/dV curve is not flat after the inelastic step and it has 
a small decay for \eV\ larger than the inelastic thresh- 
old. This non-equilibrium effect has been already ob- 
served experimentally 11 ' 27 . Using the same theory with 
a smaller tip-atom coupling (smaller vt) results in dl /dV 
lineshapes identical to those equilibrium calculations 22 . 

The non-monotonic dl/dV can be explained as fol- 
lows. As the bias goes across the inelastic threshold, 
\eV\ — 4|Z?|, there is a population transfer from the 
ground state doublet (S z — ±5/2) to the first excited 
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FIG. 2: (Color online) (a) dl/dV, in units of go, for single 
Mn in CU2N surface probed with a non-polarized tip for a 
fixed temperature ksT — 0.5K and B = 0. (b) Steady state 
populations of each eigenstate versus applied bias. Here £ = 
1, vs = 1 and vt = 1. 
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FIG. 3: Inverse of the relaxation time of the spin state S z = 
+5/2 (black dashed lines) and charging time Tq 1 = 7/jv/e 
(blue solid lines), versus the applied bias for two different 
couplings with the tip: (a) vt = 0.1, and (b) vt ~ 1. In both 
cases, X" = 0.5K, vs = 1 and f = 1. 



state doublet S z = ±3/2. This can bee seen in Fig. EJb). 
As soon as the population transfer to the first excited 
doublet takes place, a second inelastic channel opens: 
the transition from the first to the second excited state 
doublet (S z — ±1/2, whose energy is 2\D\, smaller than 
the first step). It turns out the intensity of the primary 
inelastic step (A = 4|D|), given by the matrix element 
|(±5/2|S ,± |=f3/2}| 2 , is larger than the intensity of the sec- 
ondary transition (A = 2|D|). Thus, the depletion of the 
primary transition in favor of the secondary one results 
in a decrease of the conductance. In the case of FePc 
molecules, discussed below, the secondary transition is 
stronger than the first one, resulting in an increase of the 
conductance after the first step. 

The non-equilibrium occupations can be understood as 
the balance between two driving forces. Spin-flip assisted 
tunneling events heat the atomic spin, delivering energy 
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FIG. 4: (a) Energy spectra corresponding to Hamiltonian © 
for a Mn dimer over a CU2N surface versus applied magnetic 
field. Spectrum is referred to the ground state energy and 
given in units of the exchange coupling (J = 5.9 meV). The 
magnetic field is applied in the surface plane forming a 55 
angle with the Cu-N direction. 
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FIG. 5: dl/dV curve for the Mn dimer over a CU2N surface 
proved with a non-polarized tip. Each panel corresponds to a 
pair of {v T (l),v T (2)} values: (a) {2,1.2}, (b) {0.4,0.1} and 
(d) {0.01, 0.001}. T = 0.6K, C = 1, v s (i) = 0.25. 



B. Dimer 



of the order of eV at a pace set by the inelastic current. 
The steady state is reached when the heating power is ex- 
actly compensated by dissipation. The later occurs via 
atomic spin relaxation due to exchange coupling to the 
tip and surface electrons. This process is enabled even 
at zero bias. Interestingly, the steady state occupations 
can differ enormously from the zero bias thermal equi- 
librium. At eV = 2meV, the occupation of the ground 
state doublet is half of the one in equilibrium and barely 
twice the one of the higher energy spin levels, which are 
almost empty at zero bias. 



The inverse of the lifetimes of the two competing pro- 
cesses are shown in Fig. [3] There we show the relaxation 
rate l/Ti of a magnetic spin state, i.e. S z — +5/2, as a 
function of V. When the tip is fully decoupled (no cur- 
rent through the system, 1Z = 0), the spin relaxes in a 
time scale T\ which is independent of the applied bias. 
For a small coupling, Fig. [3Ja) , the relaxation rate in- 
creases by several orders of magnitude when the bias is 
increased. This effect is even more dramatic when the 
ratio r approaches 1 [TZ = 1/2, see Fig. Hfb)]. In the 
weak coupling curve we can easily see the crossover from 
the equilibrium regime at low bias, where T\ << T q , to 
the non-equilibrium regime, for which T\ » T q . In the 
other case the crossover occurs at a much lower voltage. 
To plot these curves we take a zero bias conductance 
(7q = 0.75/iS for the vt = 1 case. 



1. Non equilibrium effects 

From the discussion above, the non-monotonic line- 
shape observed for the Mn monomer is related to non- 
equilibrium effects. Interestingly, correlation Kondo-likc 
effects could also modify the lineshape 3 -!. Given the fact 
that Kondo effect occurs in the case of a Cobalt atom 
deposited in the same surface^ this type of effect can 
not be ruled out in the Mn monomer. In contrast, the 
Mn dimer has a S = ground state d 22 ' 53 and provides 
an ideal system to test the non-equilibrium physics 2 ^. 

The Mn dimer was studied experimentally under low 
current conditions by Hirjibehedin et al£ and, more re- 
cently, under high current conditions by Loth et al^ 
They have observed a dramatic modification of the line- 
shape, which can be accounted for by our theory, as we 
show here. The Mn-Mn exchange interaction in this sys- 
tem is antiferromagnetic. The fitting- of the experimen- 
tal results to the Hamiltonian model, Eq. ([3]), gives a 
J1.2 = J = 5.9, while D, E and g are kept as for the 
monomer. 

Fig.|4]shows the lowest energy spectra of the Mn dimer. 
Since J >> \D\,E, the total spin S is a good quantum 
number at zero order in \D\/ J. Thus, the ground state 
is S = 0, the first excited state S = 1 and energy J, the 
second S = 2 and energy 3J and the third S = 3 and 
energy 6 J, all energies measured with respect to that of 
the ground state . The 25+1 degeneracy of the S > 
multiplets is weakly lifted by the small anisotropy terms 
D and E. The allowed transitions induced by the ex- 
change coupling ((3]) , when the tip is more coupled to one 
of the two atoms, satisfy AS — ±1. The lowest energy 
transitions are marked in Fig. U with vertical arrows at 
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energies Si, with i= 1,2,3. 

The experimental results of the IETS show very dif- 
ferent profiles as the current through the system is 
changed^ For low currents only the transition at energy 
81 ~ J is observed and flat plateaus appear in the dl /dV 
spectra before and after the inelastic step. This pri- 
mary step corresponds to the transition from the S = 
ground state to the first excited state 5 = 1. As the 
current is increased, by reducing the tip-atom distance, 
additional steps appear at higher energies, correspond- 
ing to the transitions between the 5=1 and 5 = 2, 
and the 5 = 2 and 5 = 3 states, with energies 82 ~ 2J 
and S3 « 3 J. In addition, the dl/dV line shapes are 
not flat away from the steps either. These results are 
reproduced by our non-equilibrium theory. Fig. [5] shows 
the theoretical dl/dV curve for three different couplings 
with the tip. When the tip is weakly coupled to the chain, 
Fig.JSJa), the step corresponding to the 5 = 0— >5=lis 
clearly visible, while excitations from the 5 = 1 — > 5 = 2 
are quenched since the 5 = 1 is only slightly populated. 
When the coupling vt is increased (higher current), tran- 
sitions 5=1^5 = 2 and 5 = 2 — > 5 = 3 become 
possible for bias \eV\ > 2 J and \eV\ > 3J respectively ( 
see Fig. El[b)). The new transitions are possible at high 
current due to a significant current induced occupation 
of the excited states 5 = 1 and 5 = 2 in the Mn dimer. 
In contrast with the Mn monomer, the the excited state 
spin flip transitions energies 82 and S3 are larger than the 
primary spin transition, resulting in in new steps in the 
spectra. These experimental results, together with the 
theoretical interpretation, provide strong evidence of the 
capability of the STM current to drive the spins of the 
magnetic adatoms. 



2. The case of symmetric coupling 

Whereas the results above are in very good agreement 
with the experimental data, 27 it is worth pointing out 
that this is only so if we assume that the exchange as- 
sisted tunneling is stronger through one of the atoms. 
However, it vanishes identically in the symmetric cou- 
pling case, wt(1) = i"r(2). In Fig. |6]we plot the height of 

/ 2 

the inelastic step, given by A s ^t = J2 a Em' S a,TS > 
as a function of the lateral position of the tip across the 
dimer axis, as modeled by the ratio wt(1)/wt(2). The 
inelastic step cancels identically when the tip is in the 
middle. This prediction of the model is at odds with un- 
published experimental data, which do not show a strong 
dependence of the inelastic current as the tip is moved 
along the Mn dimer axis. From the theoretical point of 
view, the cancellation of the exchange assisted tunneling 
in the case of the Mn dimer symmetrically coupled to the 
tip arises from the fact that, in this particular case, the 
operator in the transition matrix element (I20[) is the to- 
tal spin of the dimer, and then the eigenstates of 5 2 and 
S z are also eigenstates of V. As a result, the coupling 
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FIG. 6: Transition amplitude As^t versus the ratio 
wt(1)/«t(2). T = 0.6K, C = 1, v s {i) = 0.25. 



Hamiltonian is diagonal, and no transitions are possible. 
Notice that this problem is specific of the dimer. In the 
case of the monomer the observed spin transitions occur 
within states with the same 5 = 5/2. In the case of the 
trimer and longer chains the tip can not be coupled iden- 
tically to all the atoms and the theory accounts for the 
data^ 2 - 

There are several spin interactions other than the inter- 
atomic exchange that break the spin rotational invariance 
and could, in principle, solve the problem: the single ion 
anisotropy terms, D and E, the hyperfine coupling with 
the nuclear spin of the Mn, / = 5/2, and the direct mag- 
netic dipolar coupling. We have included them in our 
calculations, but they are much weaker than the domi- 
nant exchange, so that they do not change qualitatively 
the curve ([6]). Thus, even in spite of the apparent success 
of the perturbative approach using Eq. (U]), this partic- 
ular result indicates the presence of additional terms in 
the Hamiltonian or the need to go beyond lowest order 
in perturbation theory . Further work, going beyond the 
phenomenological theory is under way^ 



C. Magnetic molecules 

As a final example of our non-equilibrium theory 
with non-magnetic electrodes, we consider the case of 
IETS through Iron Phthalocyanine (FePc) molecules, de- 
posited on oxidized Cu surfaced FePc are flat organic 
molecules with C4 symmetry with a core made of a single 
Fe 2+ ion surrounded by 4 Nitrogen atoms embedded in 
Benzene groups. In gas phase, the crystal field of the lig- 
ands is high enough as to reduce the spin of Fe 2+ from 
5 = 2 (high spin) to 5 = 1 (intermediate spin). Be- 
cause of the C4 symmetry of the gas phase, the single 
spin Hamiltonian of the molecule has E = 0. 

According to the IETS data^ the symmetry is re- 
duced when deposited on the oxidized surface. In partic- 
ular, two adsorbed states (a and /3) were experimentally 
observed with different spin excitations^ In both cases 
the spin excitations of the FePc could be assigned to 5 = 
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FIG. 7: Differential conductance versus applied bias at zero 
magnetic field for the a (a) and j3 (b) FePc on Cu(110)(2 x 
1) — O as a function of applied bias for different magnetic field 
applied perpendicular to the sample surface, as in Ref. Il5l . 
Curves have been shifted by 0.5 units for clarity. T = 0.5K, 
Vt = 0, v s = 1, v T = 0.2 and £ = 1. 



1 but the anisotropy parameters, determined from the ex- 
perimental differential conductance curves, varied in the 
two cases. For the a (/3) configuration, D — —3.8 ± 0.04 
meV (D = -6.9 ± 0.04 meV), E = 1.0 ± 0.01 meV 
(E = 2.1±0.04meV) and g = 2.3±0.02 (g = 2.4±0.05). 
The origin of this drastic change in anisotropy deserves 
further theory work. 

The 5 = 1 single spin model can be solved analytically 
(see for instance appendix in Ref. I56I ). With D < 0, 
and E — 0, the ground state would be the S z = ±1 
doublet with energy — D below the S z = excited state. 
At finite E the ground state doublet splits, in bonding 
and anti-bonding combination of the states S z = ±1. 
The splitting is 2E. Thus, there are two spin transitions. 
The low energy one, with A = 2E and A5 Z = 0, and the 
high energy one, with energy \D\ + E and AS Z = ±1. 
The magnetic field along the z axis competes with the E 
induced splitting of the ground state. As B increases, the 
E induced mixing of the S z = ±1 components decreases, 
and so it does the primary transition, which occurs via 
A5 2 = events. 

Fig. [7] shows our non-equilibrium theoretical dl/dV 
results using the values of D and E given above. Our the- 
ory reproduces not only the evolution of the steps with 
the magnetic field, but also the mild non-equilibrium fea- 
tures reported in Ref. [TH After the first (second) step 
the conductance has a small positive (negative) slope. In 
contrast, the equilibrium theory^ yields flat steps. The 
sign of the non-equilibrium slopes depends on the relative 
value of the inelastic channel strengths of the primary 



and secondary transitions. At eV = 2E, when the first 
excited state is populated, the secondary transition, with 
energy \D\ — E becomes possible, at finite temperature. 
Since this transition between excited states has a larger 
quantum yield, the overall conductance increases. The 
opposite scenario occurs in the second step. 



IV. SPIN POLARIZED TIP 

The results of the previous section give very strong 
support to the notion that tunneling electrons can drive 
the spin of the magnetic adatoms far from equilibrium. 
Since we have been considering spin unpolarized tunnel- 
ing electrons, these non-equilibrium effects can not result 
in a net spin transfer. From the theory standpoint, this 
should change dramatically in the case of spin-polarized 
transport electrons. As it was shown in a seminal work 
by Slonczewski^ the back action of transport electrons 
on a magnetic moment can be used to rotate the magne- 
tization direction. This effect, known as spin-transfer 
torque, have been observed in nano-pillars of tens of 
nanometers^ down to tiny nanomagnets made of 100 
atoms^ but still in the semiclassical domain. In a pre- 
vious paper— we modeled the spin dynamics of a single 
Mn atom under the influence of spin polarized current. 
We found that the if the tip was spin polarized, the spin 
polarized current would result in a net spin magnetiza- 
tion of the magnetic adatom whose orientation relative to 
the tip moment would depend on the polarity of the bias, 
quite in agreement with the macroscopic spin transfer 
torque. In parallel to our work, S. Loth et. al^ demon- 
strated experimentally the single atom spin transfer. 

In our work in Ref. [25| the origin of the tip magnetiza- 
tion was ferromagnetic order. In the experiment of Loth 
et al, the spin polarized current is achieved by sticking 
a single Mn atom into the tip and applying a magnetic 
field to freeze its spin fluctuations. The external magnetic 
field affects also the surface atom. The very different role 
played by the Mn in the tip and the Mn in the surface 
underlines the important role played by spin isolation. 
Whereas it is still possible to model the spin of the Mn 
in the surface as a quantized spin weakly coupled to the 
surface electrons, this picture seems to break down for the 
case of the Mn in the tip, due to a combination of charge 
transfer, Kondo coupling and very reduced spin lifetime. 
Thus, the Mn in the tip acts as a spin filter for the trans- 
port electrons. Whereas this picture works qualitatively, 
we believe this issue deserves further worki^ 



A. Current induced spin switching 

The flow of spin polarized current through a single 
magnetic atom is expected to result in a transfer of a 
net spin into the atom. In the case of a single or a few 
magnetic atoms, where time reversal symmetry is not 
spontaneously broken at zero magnetic field , the equilib- 
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rium occupation of states with opposite S z is the same, 
resulting in a null average magnetization. Spin polar- 
ized current changes this situation via spin-flip inelastic 
tunnel^ The mechanism is the following. The dominant 
inelastic transitions in the case of the Mn monomer are: 

• Spin increasing (SI) transition, for which the Mn 
spin goes from S z = — | to S z = — § and the trans- 
port electron goes from the high energy electrode 
with spinf to the low energy electrode with spin J,. 

• Spin decreasing (SD) transition, for which the Mn 
spin goes from S z = +§ to S z = +§ and the trans- 
port electron goes from the high energy electrode 
with spinj. to the low energy electrode with spin f . 

In the case of spin unpolarized current, these two pro- 
cesses are equally likely and result in the depletion of the 
two states of the ground state doublet shown in Fig. (5Jb) . 
In the case of a spin polarized tip, the two processes are 
no longer equally likely, resulting in a net spin transfer 
from the spin current to the atomic spin. Let us con- 
sider the case where there are more 1 than t electrons 
in the tip. This means negative tip spin average (i.e., 
Vt < 0) and positive tip magnetization. When electrons 
go from the tip to the surface (V > 0) , the SD processes 
are dominant, as a result of which the positive S z states 
are depleted and a negative (S~) is expected. Thus, we 
expect that at positive bias (electrons going from tip to 
surface) the current co-polarizes the spin of the atom. 

We now consider electrons going from surface to tip. 
Since the density of states of spin 1 electrons is higher, 
the SI process is now more likely than the ST one. As 
a result, the negative S z states should be depleted, re- 
sulting in a positive atomic spin. Thus, we expect that 
V < (electrons going from surface to tip) the current 
counter-polarizes the spin of the atom. 

Our simulations confirm this scenario. We only con- 
sider the simplest case in which the tip polarization is 
assumed parallel to the Mn easy axis -perpendicular to 
the CU2N surface. We consider first the case of zero mag- 
netic field. We choose Vt < 0, because is convenient for 
the discussion at finite positive field below. In Fig. [SJa) 
we show the average atomic spin moment along the easy 
axis, as a function of the applied bias. It vanishes at 
zero bias, reflecting the absence of spontaneous time re- 
versal symmetry breaking of such a small system. At 
finite bias the magnetic moment aligns with that of the 
tip when electrons flow from tip to surface (V > 0), and 
do exactly the opposite when the electrons flow from the 
surface to the tip (V < 0). Interestingly, the average 
atomic spin is finite even when \eV\ < 4|D|, the exci- 
tation energy. This is due to the existence of thermally 
excited quasiparticles. However, the time necessary to 
drive the spin of the atom increases exponentially when 
for \eV\ < 4\D\M 

The average atomic spin increases both with the ap- 
plied voltage and the spin polarization of the tip, as given 
by Vt- The effect is null at Vt = 0, as it should, and it 
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FIG. 8: (Color online) (a) Average magnetization < S z > for 
the Mn monomer in CU2N surface probed with two differently 
polarized tips, Vt = —1/3 (red-dashed line) and Vt = — 1 
(solid-blue line) versus applied bias, (b) "Universal curve" of 
the magnetization versus tip polarization (for positive applied 
bias), (c) Steady state populations of each eigenstate versus 
applied bias. Solid lines for energy levels with S z < and 
dashed lines for S z > 0. Here £ = 0.5, vs = 1, T = 0.5K, 
B = 0. and v T = 0.7. 

is maximal for half-metallic tips Vt = ±1. For a fixed 
tip polarization the effect saturates at a certain voltage. 
Interestingly, the saturation magnetic moment depends 
only on the value of Vt and is quite independent of tem- 
perature and other parameters in the calculation. We 
discuss this universal behavior in Sec. IIV CI The non- 
zero atomic spin polarization reflects the bias induced 
breaking asymmetry of the steady state occupation of 
the two states of the ground state doublet S z = ±5/2, 
as shown in Fig.[5Jc). Notice the striking difference with 
the case of equilibrium, for which the occupations of these 
two degenerate states are identical. The steady state is 
reached thanks to the competition between bias induced 
spin-transfer and exchange induced spin relaxation dis- 
cussed in the previous section. 



B. Effects of spin polarization on transport 

Importantly, the current induced polarization of the 
atomic spin can be detected through its influence on the 
conductance of the system. The simplest effect comes 
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from the elastic magnetoresistance: conductance is larger 
when spin polarization of tip and magnetic atom are par- 
allel. In the case discussed above, this results in a larger 
conductance at large positive bias than at large negative 
bias. At small bias, there are several competing effects. 
For simplicity let us consider the case of a single mag- 
netic adatom. We can write the differential conductance 



as 



G — G\ + G2 + G3 + G4 



(26) 



where the different terms are obtained by deriving / 
in Eqs. (|22M24[) with respect to bias and are shown in 
Fig. lib): 

Gi = g (l + 2C(S z )V T ), 
G 2 = 2g V(V T d{Sz) 



G 



gs_ 



dV ' 

E p m(v) 



M.M' 



lM,M 



a 

Xi'_(A M ,w + eV) + T T Z xy (M, M')i' + (A M M> + eV) 



M.M' 



M.M' 



+P T i+{A M ,M> + eV)S xy (M,M') 



dP M (V) 
dV ' 



(27) 



with i' ± = di±/dV. E xy (M, M') is defined as in Eq. ((23)) 
but with without the weighting factors. G\ and G2 (G3 
and G4) correspond to the elastic (inelastic) contribution 
of the current. G\ gives the dominant magnetoresistive 
contribution at large bias discussed above. Gi gives a 
smaller contribution associated to the change of the av- 
erage adatom spin as a function of bias. This term is 
responsible of the non-monotonic decay of the conduc- 
tance after the pronounced change induced by G\. In 
the extreme case shown in Fig. [5J corresponding to a half 
metallic tip, G\ is the dominant contribution. Finally, 
the two inelastic contributions, G3 and G4 peak close to 
the transition energies ±4|D|. G3 corresponds to the in- 
elastic conductance, as if the occupations Pm(V) where 
bias independent, and G4 is the contribution coming from 
the fact that Pm do depend on the bias. In turn, both 
G3 and G4 have two contributions, one that is present for 
non-magnetic tip and another one proportional to the tip 
spin polarization Vt- 

Experimentally it might be hard to disentangle G2, 
G3 and G4, but not Gi which provides a direct way to 
quantify the atomic spin at large bias, denoted by ±Voc,: 



G 1 (+F 00 )-G 1 (-K ) 
G(V = 0) 



= KPriSzK+Voo) (28) 



where G(V = 0) = go is the zero bias conductance 
and we have used the fact that, at zero magnetic field, 
(S z )(+V 00 ) — —{S z )(—V 00 ). Since Gi can be the domi- 
nant contribution, replacing G\ by the total G in equa- 
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FIG. 9: (Color online) (a) Differential conductance for the Mn 
monomer versus the applied bias for two different tip polar- 
izations: Vt — — 1/3 (dashed line) and Vt = — 1 (solid line), 
(b) Each of the contributions to the dl/dV for Vt = —1/3: 
Gi (thin-black line), G2 (thick-blue line), G3 (thin-dashed 
line) and G4 (thick-green line) versus applied bias. T = 0.5K, 
£ = 0.5 and vt = 0.7. Magnetization direction was fix parallel 
to the easy axis and no magnetic field was applied. 
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FIG. 10: (Color online) Average magnetization {&) for the 
Mn monomer versus applied bias for B — 7T corresponding 
to the low current (blue line) and flight current (black line) 
regimes. The tip magnetization direction was fix parallel to 
the easy axis and to the applied field. Vt = — 0.31,T = 0.5K 
and vs = 1. 



tion (|28[) can give a rough estimate of the quantities in 
the right hand side of that equation. 



1. Finite magnetic field 

We now analyze how the magnetic field used in the ex- 
perimental setup^l changes the picture discussed above. 
In these experiments, the tip polarization is achieved by 
attaching a single Mn atom at the tip apex and applying 
a very intense magnetic field (7T) perpendicular to the 
surface. In the experiments with spin polarized tips, 27 
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FIG. 11: (Color online) (a) Differential conductance for B = 
7T versus applied bias for the low current regime, (b) Each 
of the contributions to the dl/dV: Gj (thin-black line), Gi 
(thick-blue line), G3 (thin-dashed line) and G2 (thick-green 
line) versus applied bias. T = 0.5K, £ = 0.5 and vt = 0.7. 
More detail in the text. 
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FIG. 12: (Color online) (a) Differential conductance for B = 
7T versus applied bias for the hight current regime, (b) Each 
of the contributions to the dl/dV: Gi (thin-black line), G2 
(thick- blue line), G3 (thin-dashed line) and G2 (thick-green 
line) versus applied bias. T = 0.5K, £ = 0.5 and vt = 0.7. 
More detail in the text. 



the dl/dV curve changes radically from low current to 
high current . In both cases, the dl/dV curves are not 
even with respect to bias, as expected from the discus- 
sion in the previous section. However, the lineshapes 
differ significantly from the ones described in our pre- 
vious work2£ and in the previous section (Fig. The 
origin of this discrepancy can be traced back to the ef- 
fect of the magnetic field on the surface adatom, absent 
in our previous calculation. Since k^T « g/i^B, the ap- 
plied field already polarizes completely the atomic spin 
at zero bias. In particular, this means that the occupa- 
tion of the atomic spin state S z = —5/2 is very close to 1 
and at zero bias the average spin of the adatom is finite, 
negative, and parallel to the tip spin polarization, which 
maximizes the G\ term in the conductance. 

When bias is applied, the occupation of the ground 
state S z — —5/2 is depleted, which degrades the elastic 
contributions to the conductance. This is reflected in the 
evolution of the average atomic spin as a function of bias, 
shown in Fig. [TU] for two cases: low current (vt = 0.08, 
C = 0.44) and high current(i> T = 0.7, C = 0.5). The val- 
ues of vt and £, as well as the value of Vt = 0.31, where 
chosen to reproduce the experimental conductance data 
of Loth et alwL. It is apparent that in both cases the de- 
pletion of the average spin is larger for negative bias, for 
which transport electrons tend to counter-polarize the 
tip, than for positive bias, for which the depletion can 
be interpreted as non-equilibrium heating of the atomic 
spin. In the large current case, the current is able to re- 
verse the average atomic spin from the equilibrium —2.5 
at zero bias to +1 at -20 mV. The corresponding change 
in the elastic components of conductance, G\ and G2, due 
to spin contrast, is shown in figures (jllb ) for low current 
and in (|12b ) for high current. In the high current case, 
these contribution are dominant, and explain the exper- 
imental observation that the conductance is smaller at 
negative voltage. Thus, in the high current case the fact 
that G(+20) < G(— 20) can be linked to the reversal of 
the atomic spin at negative bias and the resulting reduc- 



tion of the conductance, compared to zero bias. 

In the low current case the small changes in the atomic 
spin make G\ and G2, the elastic magnetoresistive con- 
tributions, quite independent of bias, so that change as a 
function bias so that they simply offset the total conduc- 
tance. This results in a conductance which is minimal 
at zero bias and the larger at large negative bias than 
at large positive bias. As opposed to the high current 
case, where the asymmetry is dominated by Gi and G2, 
reflecting the reversal of the average atomic spin, in the 
low current case the asymmetry comes mostly from the 
asymmetry in the inelastic step G3, associated to the 
term proportional to Vt- 



C. Universal magnetization profile 

In Fig. |5fc) we plot the saturation magnetization, 
reached after a large enough bias is applied, as a function 
of tip polarization. We have verified that this curve turns 
out to be independent of all transport parameters, C, r, 
v n ... and it depends only on the spin of the magnetic 
atom S and anisotropy parameters D and E. In fact, it 
does not depends either on the temperature, as long as 
our approach of neglecting the phonon contribution re- 
mains valid. When combined with eq. (|2"8")l . this could 
be used to determine the tip polarization. 

The universality of the saturation atomic magnetiza- 
tion comes from the fact that, at large bias, z±(Am,m' + 
eV) ~ ±eV so that the equation (|A3|) for the rates can 
be simplified to: 



W M 



M' 



oc V 



?M,M'\2 



Pt 



PTl igM,M'|2 

Pt 



(29) 



Making use of Eq. (|2"9"|) . the master equation for the 
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occupation of the spin states in steady state reads 



The main conclusions are the following: 



E 



M' 



\s 



M,M 1 2 



1 



^£| 5 M,M< |2 



2 

x {P M - - Pm) = 0, V M, (30) 

where the prime indicates sum is done over M' ^ M . 
Eq. (f30f shows that, in the large bias limit, the atomic 
spin steady state occupations Pm, and consequently the 
average magnetization {S z ), depend only on the matrix 
elements of the spin operators, and the polarization of 
the tip, and do not depend on the coupling strength to 
the tip and surface. 

V. SUMMARY AND CONCLUSIONS 

We have studied the mutual influence of non- 
equilibrium transport electrons and the spin of one and 
two magnetic adatoms in STM configuration. Our re- 
sults indicate that non-equilibrium effects are essen- 
tial to understand present IETS STM spectra of mag- 
netic adatoms. Our theory is able to describe correctly 
the experimental observations of IETS on single Mn 
atoms, both with non-magnetioi^and magnetic^ tips, 
on Mn dimers at low and high current^!, and on FePc 
molecules^. Our theory is based on a phenomenological 
spin model to describe both the atomic spin states and 
their coupling to the transport electrons. Current is cal- 
culated to lowest order in the tunneling Hamiltonian and 
depends on the occupation of the spin states Pm- The 
spin dynamics is described with a master equation for 
the Pm that includes the effect of carrier induced spin 
relaxation and the spin pumping due inelastic spin-flip 
assisted tunneling. 

When the time elapsed between inelastic spin flip 
events, T q is much longer than the atomic spin relax- 
ation time, Ti, the atomic spin remains in equilibrium 
and Pm does not depend on bias. In that limit, the 
differential conductance dl/dV is piecewise constant 22 
except for the steps when the bias voltage matches the 
energy of the spin excitations. When the T q is compara- 
ble or smaller than T\ , the occupations of the atomic spin 
states are driven away from equilibrium and they depend 
on the bias voltage. Whereas non-monotonic dl/dV had 
been already observed experimentally, the recent results 
reported by Loth et al£Z- have confirmed this scenario by 
controlling the tip-adatom distance. In addition, the use 
of spin-polarized tips amplifies the changes in the dl/dV 
curves as the conductance is increased. 

The results of Loth et al£Z- also indicate that, in the 
case of magnetic tips, the orientation of the average 
atomic spin can be switched at will from parallel (V > 0) 
to antiparallel (V < 0) with respect to the magnetic 
tip. The control of the spin of a single atom and a sin- 
gle magnetic molecule had been predicted by theor y 22 i 46 
The control of atomic spin with non-equilibrium spin 
polarized carriers is similar to that obtained by optical 
pumping^. 



1. The dynamics of the atomic spin under the influ- 
ence of tunneling electrons is governed by two in- 
trinsic time scales, the inelastic transport time T q 
and the atomic spin relaxation time T\. When the 
current induced spin flips occur more often than the 
time it takes to the atomic spin to relax, i.e., when 
T q < Ti, non-equilibrium effects build-up. This 
makes the occupation of the spin states different 
from that of equilibrium. 

2. The conductance lineshape is sensitive to the oc- 
cupation of the atomic spin states. This might be 
used to perform transport-detected single atom res- 
onance experiments. Our calculations indicate that 
non-equilibrium effects have been observed in single 
Mn atoms, in Mn dimers and in FePc molecules. 

3. A rough estimate of the quality factor of the spin 
excitation with energy A = Huj, defined as Q = 
luTi, can be obtained from the transport experi- 
ments. We assume that the current at which the 
low temperature conductance line shape starts to 
deviate from a piecewise constant function yields 
T q ~ T\. Let AGi n and V be the height and bias 
of the primary inelastic step. Then, the inelastic 
current is /in = AGqV — e/T q . Then we get 



ujT\ 



AG' 



(31) 



4. Spin polarized STM can be used to magnetize the 
atomic spin both parallel or antiparallel to the mag- 
netic tip moment. When electrons tunnel through 
the magnetic atom from the magnetic tip to the 
surface, the atomic is magnetized parallel to the 
tip. Reversing the bias results in an opposite spin 
polarization. 

5. The bias induced adatom spin polarization results 
in asymmetric conductance lineshapes due, in most 
part, to the dependence of conductance on the rel- 
ative orientation of the adatom and tip magnetiza- 
tions. 

6. The saturation atomic spin magnetization, ob- 
tained at large bias, is only a function of the impu- 
rity spin, the anisotropy parameters and the polar- 
ization of the tip. 

Future work should address open problems, like the 
origin of the spin assisted tunneling Hamiltonian for spin 
larger than 1/2, the effect of atomic spin coherence, which 
we have neglected in the master equation ((5]) and the fact 
that the observed inelastic steps in the Mn dimer do not 
depend on the lateral position of the tip, in contrast with 
our theory. 
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Appendix A: Equation for the rates 

In this appendix we derive the general expression of the 
transition rates. Applying the Fermi Golden rule using 
the tunneling Hamiltonian (j4|) as perturbation, one gets 



-k+(W) 



iM.Af' 



(A4) 



Here we have introduced TZ ± (rjrj l ) = p v ^p V '-t ± Pr^/Vi 
and the operators S± = S x ±iS y . Notice that expression 
(|A4j) is valid for all states M, M' , in contrast to Eq. (5) 
of Ref. 25 where only the inelastic matrix elements were 
explicitly written. 



Appendix B: Equations for the current 



1 kaM.k'a'M 1 



2n 



S (e ( T7 7 (fc) + E M - e a ' V >{k') - E M >) 



,(«) 



Erp ' era' qM, 
±a 9 °a,r, 



M 



(Al) 



The modulus square in Eq. (| Al|) can be expanded to 
obtain 



E T ^ 



T (a ) T iP) 



.t,t ■ ,T'g gM,M' gAf',Af 



a, 777/ ' 



a .p .(7 ct 



Considering the explicit form of the Pauli matrix ele- 
ments, the sum over a and a 1 can be done, with just a 
few non-zero contributions. Using the definition of the 
transition rates W^J,, Eq. ([7]), we obtain after some 
algebra 



Here we shall derive expressions (fT5)l and (|T5)) . Let 
us start with the expression for the current, Eq. (|16p . 
Dividing the contribution into its elastic and inelastic 
part, with the help of Eq. (|A3[) and (|A4I) . we can write 



LMR 



^eT oXP s PT [1 + 2C{s ^ s)VT] 

x (t(eV) + G{-eV)) , (Bl) 



(A2) or i using the definition of go and i+ 



Iq + Imr = -[l + 2{(S ZtT>s )V T ]i-(eV) ) 
e 



which is the result of Eq. (1181) . For the inelastic contri- 
bution (M t^M'), the difference - W^fXf, can be 
written with the help of Eqs. (|A3|) and (IA4|) as 



vv M,M' 



h 



-g(A M ,M' + H - Mn') E M,jUA3) 
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Using the definitions of g s and S xy (M, M'), expression 
(fl9t is recovered. 
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